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ABSTRACT 

The 2-Micron All Sky Survey (2MASS) (Skrutskie et al. 1997) and the DEep Near 
Infrared Survey of the southern sky (DENIS) (Epchtein et al. 1994) have revealed a 
heretofore unknown population of free brown dwarfs that has extended the local mass 
function down to as small as O.OlMg (Reid et al. 1999). If this local proportion of brown 
dwarfs extends throughout the Galaxy — in particular in the Galactic bulge — one expects 
an increase in the predicted fraction of short time scale microlensing events in directions 
toward the Galactic bulge. Zhao et al. (1996) have indicated that a mass function with 
30-60% of the lens mass in brown dwarfs is not consistent with empirical microlensing 
data. Here we show that even the much lower mass fraction (~ 10%) of brown dwarfs 
inferred from the new discoveries appears inconsistent with the data. The added brown 
dwarfs do indeed increase the expected number of short time scale events, but they 
appear to drive the peak in the time scale frequency distribution to time scales smaller 
than that observed, and do not otherwise match the observed distribution. A reasonably 
good match to the empirical data (Alcock et al. 1996) is obtained by increasing the 
fraction of stars in the range 0.08 < m < 0.7Mq considerably above that deduced from 
several star counts. However, all inferences from microlensing about the appropriate 
stellar mass function must be qualified by the meagerness of the microlensing data and 
the uncertainties in the Galactic model. 



Subject headings: Galaxy: stellar content — stars: low-mass, brown dwarfs 



1. Definitions and Assumptions 

A microlensing event is said to occur whenever a 
source star and lens star pass each other at an angu- 
lar separation that places the source within an Ein- 
stein ring radius of the lens. The gravitational fo- 
cusing of the source light by the lens leads to the 
observed amplification of the source brightness. The 
time scale of the event is defined as — Re/v, where 
Re = ^AGmDor.{Dos - Doe)/{c^Dos) is the Ein- 
stein Ring radius with Dql and Dqs being observer- 
lens and observer-source distances respectively, m the 
mass of the lens, G the gravitational constant, c the 
velocity of light, and v is the relative transverse ve- 
locity between lens and source projected onto the lens 
plane. Empirically, the time scale is approximately 
half of the duration of an event. 

To illustrate the effect of a population of brown 
dwarfs on the microlensing time scale frequency dis- 
tribution, we adopt the model of Zhao (1996) for the 
mass distribution in a bar-like Galactic bulge. 
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with (7„ — 0.6, cr^ = 1.49 kpc, cfy = 0.58 kpc and 
= 0.40 kpc. Here x,y,z' are orthogonal coor- 
dinates along the principal axes of the triaxial bar. 
The coefficient p^„ is chosen such that the mass of the 
bulge is 2.2 x IO^^Mq (Zhao 1996). To this we add the 
exponential disk model of Bahcall & Soneira (1980) 
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where z' is the coordinate perpendicular to the plane 
of the Galaxy and r is the radial coordinate in the 
plane of the Galaxy. We choose — 2.7 kpc for 
the scale length in the radial direction (Zhao, Spergel 
& Rich 1995; Kent, Dame & Fasio 1995) instead of 
3.5 kpc chosen by Bahcall & Soneira. The coefficient 
Pao is chosen such that = 0.05Mq/pc'^ at z' = 
and r = 8 kpc, the Sun's distance from the Galactic 
center. The observations are consistent with the long 
axis of the bulge inclined about 13° to 20° relative to 



the line of sight to the Galactic center with the near 
side of the bar lying in the first quadrant (Zhao 1998). 
An inclination of 13° is chosen in the following. 

We normalize all distances by Ds — 8 kpc, all ve- 
locity components by Wlsr — 210 km/sec, the circu- 
lar velocity of the local standard of rest (LSR), the 
masses by the solar mass Mg and the timescale 
by to = V'iGMoiV^SIc^ = 66.72 days, which 
is the timescale for an event with a solar mass lens 
located at a distance of 4 kpc with a source at rest 
at the Galactic center (8kpc). For to, the lens ve- 
locity and observer velocity are both assumed to be 
the circular velocity Wlsr — 210 km/sec in the az- 
imuthal direction consistent with a fiat Galaxy rota- 
tion curve. With these normalizations, the number 
of events s~^ source^^ (unitt^;)^^ at timescale ts is 
(Peale 1998) 
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where <e appears only in the expression for m in 
dpL (z, m) I dm with 
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In Eqs. (^ and (||) v^, are the components of 
the relative transverse velocity v in the directions 
of increasing Galactic latitude and longitude respec- 
tively, pL{z,m) is the lens mass density, where the 
numerical coefficients in Eq. ^ mean it is ex- 
pressed in Mg/pc^, ( = Dos/Ds, z = DolIDs,. 
The distributions of random velocities about the cir- 
cular velocity of the Galactic model are /^^ = 

exp (-i;fy(2o-JJ)/(^27rcrx„), where axy is the rms 
value of the random velocity v'^^ and the subscript X 
is either L ov S for lens or source and the subscript y 
is 6 or £ for the latitude and longitude directions. The 
spatial density of sources is [C) , and the exponent 
2 -f 2/3 follows from the increase in area of the cross 
section of the angular field with distance and from 
the assumption that the fraction of stars with lumi- 
nosities greater than some L» varies as (Kiraga & 
Paczyhski 1994). 
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The details of the derivation of Eq. (||) are given 
in Peale (1998). The expression for the rate of events 
for a single source at Dqs that has random veloc- 
ity components within dv'gf^dv'^^ of {v'sf,,v'g^) for lenses 
within dDoL of Dql with random velocities within 
dv'^^dv'^i of {v'^^,v'^^), and with lens mass within dm 
of m is multiplied by the fraction of sources within 
dDos of Dqs and by 5{tE — t'^) to select a particu- 
lar time scale. Independent variable m is replaced by 
tE and the expression integrated over the ranges of all 
the variables as a function of t'^ to yield the time scale 
frequency distribution of Eq. (^) . The integral limits 
correspond to the variables in the same order as the 
differentials in the integrand, where the integration 
over C starts at 0.1 instead of zero to avoid a singular- 
ity. This is justified as there are not likely to be any 
sources involved in microlensing events closer to us 
than 0.8 kpc. The distributions of lenses and sources 
along a particular line of sight are identical and fol- 
low from the above Galactic model. The assumed cir- 
cular velocities for the general Galactic rotation are 
Vlo, Vso = 210 km/sec for z, C < 0.5875 (r > 3.3 kpc) 
consistent with a flat rotation curve except the circu- 
lar velocities decrease linearly to at the Galactic cen- 
ter for 0.5875 < z, C < 1-0 (r < 3.3 kpc). The velocity 
dispersions in the disk are {aat, (Jde) = (16, 20) km/sec 
{z, C < 0.5875), whereas the bulge velocity dispersions 
are (a,,, a,,) = (110,110) km/sec (0.5875 < z,C < 
1.0). The velocity dispersions <Ti,i,aLi, asb, <^se are as- 
signed to disk or bulge dispersions in Eq. (Q) accord- 
ing to the locations of the lens and source. Evalua- 
tion of the integral with a Monte Carlo technique is 
described in Peale (1998). 

2. Results 

Estimating the mass function for brown dwarfs is 
a complicated process, since, unlike main sequence 
stars, there is no unique correspondence between mass 
and spectral type. The spectral type L dwarfs, obser- 
vations of which are the basis of the Reid et al. (1999) 
discussion, have effective temperatures between 1400 
and 2000 K with no CH4 absorption, but every brown 
dwarf with 0.01 < m < 0.075Mq passes through the 
L dwarf region as it cools, with those of higher mass 
passing through later (Burrows, et al. 1997). The 
number of L-type brown dwarfs at any epoch depends 
on the birth rate, the initial mass function, and the 
evolution of these failed stars (Reid, et al. 1999). If 
the birth rate is assumed known and a theory of evo- 
lution adopted, the dependence of the number of visi- 



ble L dwarfs on the mass function can be used to con- 
strain the latter. If a power law distribution of mass is 
assumed, the steeper distributions lead to more visible 
L dwarfs, and the index is constrained by the observed 
volume density of L dwarfs in the region reasonably 
close (< 8 pc) to the Sun. This procedure leads to 
dn/dm ^ m"^ '^ for 0.01 < m < 0.08Mq based on in- 
formation from the DENIS and 2MASS surveys (Reid 
et al. 1999). Although model dependent, this result 
is the first determination of the brown dwarf mass 
function based on observational data. What are the 
implications for the microlensing time scale frequency 
distribution? 

Several star counts yield dn/dm ~ m~^ for 0.08 < 
TO <~ 0.7— I.OMq with steeper slopes for more mas- 
sive stars (e.g. Basu & Rana 1992; Holtzman et al. 
1998; Reid et al. 1999). We shall use the Holtzman 
et al. mass functions and modifications thereof in the 
examples. Figure 1 shows the time scale frequency 
distributions for the above Galactic model for a line 
of sight toward {t,h) = (1°,— 4°), which is represen- 
tative of an average (see Peale 1998) over the direc- 
tions to 24 fields observed by the MACHO (MAssive 
Compact Halo Objects) group during the 1993 bulge 
season. During this season 12.6 million stars were 
monitored for 189 days (Alcock et al. 1996) leading 
to the choice of ordinate in Figure 1. The distribu- 
tions are shown for the Holtzman et al. (1998) mass 
function (solid line) and three modifications, two of 
which include brown dwarfs down to O.OlMg to ac- 
commodate the new discoveries. For continuation of 
the main sequence mass function into the brown dwarf 
region with slope -1, the mass of the brown dwarfs is 
only 5.5% of that of the main sequence stars, and 
it is about 7% of the main sequence star mass with 
the Reid et al. slope of -1.3, where continuity of the 
mass function is assumed. These brown dwarf masses 
are 8% and 10%, respectively, of those main sequence 
stars with 0.08 < to < I.OMq — a range more appro- 
priate to the solar neighborhood. The third modifi- 
cation contains no brown dwarfs but has many more 
low mass, hydrogen burning stars, and it yields a time 
scale frequency distribution closely matching the em- 
pirical data. (The peaks of these curves are raised by 
about 10% if the value oIvi^sr = 270 km/sec found by 
Mendez et al. (1999) is appropriate for the flat part 
of the Galaxy rotation curve.) All the mass functions 
are continuous and are normalized to yield a bulge 
mass of 2.2 x lO^'^Mg and a mass density at the Sun's 
position in the disk of 0.05Mq/pc^. The empirical 
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time scale frequency distribution from the MACHO 
1993 bulge season (Alcock et al. 1996), adjusted for 
observing efBciencies and binned into 5 day intervals, 
is shown for comparison. 



. 1 ' ' 


' ' 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 

(/,b)=(r,-4-) 


1 i 




dtg=Number of events 
with tj,,<t;,<tj.2 for 
12.6 million stars. 


- i ' 

' I' 
"i 

ij 




m->, 0.08<m<0.7 

m-2=, 0.7<m<2.0 

m-', 0.01<m<0.7 

0.7<m<2.0 

m-i-3, 0.01<m<0.08 

m-', 0.08<m<0.7 
m-2=, 0.7<m<2.0 

m-=^ 0.08<m<2.0 

MACHO 93blg 




. . 1 


. . . , \ 



0.5 1 1.5 2 

tj, (66.72 days) 



Fig. 1. — Comparison of predicted microlensing time 
scale frequency distributions with the MACHO ob- 
servational data corrected for observing efBciencies. 
The Zhao (1996) density distribution of the triaxial 
galactic bulge and the Bahcall-Soneira (1980) density 
distribution of the disk (with disk scale length 2.7 
kpc) are assumed. The average rotation curve of the 
galaxy is assumed flat at 210 km/sec for galactocen- 
tric distance r > 3.3 kpc and declining linearly to 
zero at the galactic center for r < 3.3 kpc. Veloc- 
ity dispersions in the directions of increasing galac- 
tic latitude and longitude about the circular motion 
are {(7^^,(7^,) = (16,20) km/sec for r > 3.3 kpc and 
{(^bb, f (,<) = (110, 110) km/sec for r < 3.3 kpc. 

As found by several authors (e.g. Peale 1998), mass 

functions with dn/dm ^ in the M star region 

fail to yield a sufficient number of short time scale 
events near the peak of the empirical distribution for 
the simplest Galactic models consistent with observa- 
tions. Adding the Reid et al. brown dwarfs to the 
Holtzman et al. distribution with either the -1.3 or -1 
slope increases the number of short time scale events 
as expected, but both extensions of the mass function 
fail to yield a sufficient number of short time scale 
events near the peak of the Alcock et al. results, and. 



perhaps more importantly, the peaks of the predicted 
distributions are shifted toward smaller time scales (3- 
5 days) than those of the observed peak (5-10 days). 
(See Peale (1998) for similar consequences for other 
extensions of the mass function into the brown dwarf 
region.) A fairly good match to the empirical data is 
obtained in Figure 1 by increasing the fraction of low 
mass, but still hydrogen burning stars with a contin- 
uation of the -2.2 slope of the Holtzman et al. mass 
function down to m = 0.08Mq. 

3. Discussion 

The chosen model of the Galaxy used in calculat- 
ing the predicted time scale frequency distributions 
is consistent with observations (Zhao 1996, 1998), 
although many parameters are poorly constrained 
(Zhao 1998). If the Galaxy model is assumed to be 
close to the real Galaxy, and the local population of 
brown dwarfs extends throughout the Galaxy, then 
the MACHO collaboration should have seen a larger 
fraction of events with 3-5 day time scales. One could 
conclude that the absence of these events means that 
there are fewer brown dwarfs in the bulge, where; most 
of the lenses reside, than observed locally. However, 
the same caveats and reservations about microlensing 
constraints on the mass function emphasized by Peale 
(1998) apply here. 

The meagerness of the published microlensing data 
for bulge events precludes definitive microlensing con- 
straints on the distribution of stars in the Galaxy or 
the mass function, although there are more than 200 
events that have not yet been analyzed that would 
improve this situation. We also point out that the 10 
events in the bin corresponding to time scales from 
5 to 10 days in Figure 1 involve two clump giant 
sources with detection efficiencies of about 0.6 and 
eight much less luminous turnoff and possibly more 
blended stars as sources with detection efficiencies of 
0.15 to 0.2 (Alcock et al. 1996). The average of the 
reciprocals of the efficiencies is nearly 5 meaning this 
element of the histogram is almost 5 times higher than 
the actual number of detected events. The sampling 
efficiencies for a given time scale event for a distribu- 
tion of source brightnesses are determined empirically 
and should be reasonably sound, but the reduction of 
these efficiencies due to blending of the source stars 
with other nearby stars by a factor of 0.75 is just 
an estimate (Alcock et al 1996). Until the detection 
efficiencies for bulge events are determined more rig- 
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orously, even the empirieal distribution of MACHO 
time scales that we are trying to match is only a first 
approximation. Could the detection efficiency for the 
3-5 day events be even smaller than the Alcock et al. 
estimate? If so, this element of the histogram would 
be raised, although it seems unlikely that this region 
could become the peak of the distribution and thereby 
accommodate the brown dwarfs. 

If the empirical distribution approximates reality, 
the brown dwarfs might be accommodated by reduc- 
ing the velocity dispersion (See examples in Peale 
1998), but unrealistically low dispersions would be 
necessary to shift the peak to the MACHO value, 
and increasing the mass of the Galaxy to yield the 
observed peak frequency of events would result in too 
many events with ts ^ 20 days. Finally, one notices 
the several very long time scale events in the MA- 
CHO distribution that are not consistent with any of 
the models. These events may be due to there being 
more close lenses, such as a spiral arm concentration, 
where the relative angular velocity of source and lens 
is small, than are accounted for in the axisymmet- 
ric disk model. Other than the triaxial bulge, devi- 
ations from axial symmetry in the Galaxy model are 
ignored. In addition, streaming velocities superposed 
on the Galactic rotation are not accounted for. Could 
some combination of such additional degrees of free- 
dom hide a Galaxy- wide distribution of brown dwarfs 
in the microlensing data? One cannot be sure with- 
out better constraints on Galactic structure and ve- 
locity parameters. For the present, the addition of the 
brown dwarfs to the mass function appears to widen 
the discrepancies between predicted and observed mi- 
crolensing time scale frequency distributions. 
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